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Summary. We consider a one-dimensional stochastic reaction-diffusion generaliz- 
ing the totally asymmetric simple exclusion process, and aiming at describing single 
^ lane roads with vehicles that can change speed. To each particle is associated a 

C"j jump rate, and the particular dynamics that we choose (based on 3-sites patterns) 

ensures that clusters of occupied sites are of uniform jump rate. When this model 
is set on a circle or an infinite line, classical arguments allow to map it to a linear 
network of queues (a zero-range process in theoretical physics parlance) with expe- 
ct nential service times, but with a twist: the service rate remains constant during a 

busy period, but can change at renewal events. We use the tools of queueing the- 
ory to compute the fundamental diagram of the traffic, and show the effects of a 
condensation mechanism. 
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1 A multi-speed exclusion process 

The totally asymmetric exclusion process (tasep) is a popular statistical 
physics model of one-dimensional interacting particles particularly adapted 
to traffic modeling. This is due to its simple definition, and to the non-trivial 
exact solutions which have been unveiled in the stationary regime [1]. One im- 
portant shortcoming of this model is that it does not allow particles to move 
t-H at different speeds. Cellular automata like the Nagel-Schreckenberg model [2] 

address this issue, leading to very realistic though still simple simulators. How- 
ever, these models are difficult to handle mathematically beyond the mean 
field approximation [3] and an approximate mapping with the asymmetric 
chipping model suggests that the jamming phenomenon takes place as a broad 
^ crossover rather than a phase transition [4]. In this paper, we are interested 

in analyzing the nature of fluctuations in the fundamental diagram (fd), that 
is the mean flow of vehicles plotted against the traffic density. To address this 
question, we propose to extend the TASEP in a different way, more convenient 
for the analysis albeit less realistic from the point of view of traffic. 
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Fig. 1. The basic 1-dimensional tasep model. 
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Fig. 2. The multi-speed 1-dimensional tasep model. 



The elementary 1-dimensional tasep (Fig. 1) model is defined on a discrete 
lattice (e.g. a finite ring with boundary periodic conditions, a segment with 
edges or an infinite line), where each site may be occupied with at most one 
particle. Each particle moves independently to the next site (say, to the right), 
at the times of a Poisson process with intensity \i. Therefore, the model is a 
continuous time Markov process, which state is the binary encoded sequence 
at € {0,V} N of size AT (the size of the system), where the letter V (resp. 
O) denotes a vehicle (resp. an empty space) at site i. Each transition involves 
two consecutive letters when a particle moves from site i to site i + V. 

VO 4 OV. 

In order to encode various speed levels, we propose to extend the basic 
tasep by allowing the particle to jump at different possible rates which them- 
selves vary stochastically in time (Fig. 2). Assuming for now a finite number 
of n — 1 speed levels, the Markov chain that we consider is a sequence encoded 
into a n-alphabet {O, A,B,.. .} 

a t = {Vi,...,i = l...N}, Vi€{0,A,B,...} 

where O is again an empty site, and A, B, . . . represent occupied sites with 
jump rates /x a ,/i(,, .... In our model, the transitions remain local: the particle 
may jump to the next site only if it is empty, we allow the final state to 
be conditioned by the site after the next. More precisely, we assume that 
any transition involves three consecutive letters, and distinguish between two 
cases: 

...XOY ..M...OYY Y^O rule 1, 

...XOO ...^ ...OZO Z^O rule 2. 

In the second case, the type (or equivalently the jump rate) of the particle 
is chosen randomly according to a distribution F. As a limiting case, we will 
consider a general continuous distribution F on R + . In other words, a particle 
at site i with rate [i x jumps to site i + 1 and acquires a new rate fi z which is 
a random function of Vi+2. 
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The basic assumption is that if a car gets in close contact to another 
one, it will adopt its rate. Conversely, if it arrives at a site not in contact 
with any other car, the new rate will be freely determined according to some 
random distribution. This models the acceleration or deceleration process in 
an admittedly crude manner. This setting is different from usual exclusion 
processes with multi-type particles, each having its own jump rate. It is more 
in line with the Nagel-Schreckenberg model, with the difference that only local 
jumps are allowed and speed is replaced by jump rate. 

2 X-stage tandem queue reformulation 

In the context of exclusion processes, jams are represented as cluster of par- 
ticles. Clustering phenomena can be analyzed in some cases by mapping the 
process to a tandem queueing network (i.e. a zero range processes in statis- 
tical physics terms). For the simple tasep on a ring two dual mappings are 
possible: 

• the queues are associated to empty sites and the clients are the particles 
in contact behind this site, 

• the queues are associated with particles and the clients are the empty sites 
in front of this particle. 

By using one of these mappings, the tasep is equivalent to a closed cyclic 
queueing network, with fixed service rates equal to the jumping rate [i of 
the particles. Steady states of such queueing network have been analyzed 
thoroughly (see for ex. Kelly [5]) in terms of a simple product form structure 
which we expose now. 

Consider an open L-stage tandem queue, with arrival rate A and a common 
service rate [i: L queues with service rate \i are arranged in successive order 
(the departures from a given queue coincide with the arrivals to the next one) 
and the arrival process of the first queue is Poisson with intensity A. Each 
queue is stable when A < /j,, transient when A > fi. It is then well known that 
the distribution of the number of clients X\, . . . , Xl in the queues is 

L 

P({X i =x i ,i = l...L}) = l[P x (X i =x i ), (1) 

where 

P x (X i = x) = (l-p)p x , with p d = f -. 

A 4 

If the network is closed (the last queue is connected to the first one in the 
ring geometry), then expression (1) remains valid, with the constraint that 
the total number of clients ^2 i=1 Xi is fixed. In this case, A can be chosen 
arbitrarily, as long as each queue in isolation remains ergodic. 
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Fig. 3. Mapping of the variable-speed tasep to a tandem queue 

It can be shown easily [6] that, for a plain tasep on a ring, the size of the 
jams is asymptotically a geometric random variable with parameter p (when 
N = kL. k = 1,2...). In the open geometry, the arrival rate is an external 
parameter which can be set between and p. When it becomes comparable 
to the service rate, i.e. when p ~ 1, large queues may form and a random walk 
first time return calculation yields a realistic scaling behavior for the lifetime 
distribution of jams [7] 

P(t)~t- 3/2 . 



In our multi-speed exclusion process, particles are guaranteed by construc- 
tion to form clusters with homogeneous speed, and the mapping of empty sites 
to queues is suitable (Fig. 3). The new feature is that the service rate Ri of 
a given queue can change with time: it is drawn randomly from a distribu- 
tion with cumulative distribution F when the first customer arrives. It is 
assumed that there exists a minimal service rate /io > such that F(po) = 0. 
po = X/ po < 1 is therefore the maximal possible load. The state is determined 
by the pair (X, R) and the possible transitions are as follows: 
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Since these transitions form a tree (see Fig. 4(a)), each queue in isolation 
is a reversible Markov process and its stationary distribution reads: 

P X {X = x,Redp) = P X (X = O)F(dp) Q j , 

with 
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Fig. 4. State graph for isolated queues in the case of (a) the multi-speed process 
definition, and (b) a non-reversible queue with hysteresis. States are represented by 
black dots and transitions by arrows. 



The distribution of the number of customers in the queue is therefore no 
longer geometric: 



Px(X = x) 



P X (X = x,Red^) = P X (X = 0) / F{dfi) - 

Mo J Ho VM 



(2) 



Nevertheless, the product form expression (1) for the invariant measure 
remains valid, because of the reversibility of the individual queues taken in 
isolation (see again [5]). The stationary distribution of the L-stage tandem 
queue takes the form 

L 

Px(S) = l[Px(X i ,R i ), 

i=l 

for any sequence S — {(Xi,Ri),i = 1, . . . , L}. While each queue has a different 
service rate at a given time, all the queues have globally the same distribution. 
Our model is therefore encoded in the single queue stationary distribution 
Px(Xi,Ri). 



3 The fundamental diagram 

As announced in Section 1, we turn now to the fundamental diagram (fd), that 
is the plot of the mean flow of vehicles against the traffic density. By nature, 
the fluctuations in the FD are associated to the jam formation. Schematically, 
three main distinct regimes or traffic phases have been identified by empirical 
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studies [8]: one for to free-flow, and two congested states, the "synchronized 
flow" and the "wide moving jam". 

In the case of the basic tasep, it is well known, and rigorously proved in 
some cases, that an hydrodynamic limit can be obtained by rescaling both 
the spatial variable x = i/N and the jumping rate according to p(N) — 
NVo, where N is a rescaling which we let to oo and Vb is a constant. The 
corresponding coarse grained density p satisfies the inviscid Burger equation 

The FD at this scale is deterministic, since 

J( P ) = V oP (l - p), 

and symmetric w.r.t. p = 1/2 because of the particle-hole symmetry. Vb rep- 
resents the free velocity of cars, when the density is very low. 

In practice, points plotted in experimental fd studies are obtained by 
averaging data from static loop detectors over a few minutes (see e.g. [8]). This 
is difficult to do with our queue-based model, for which a space average is much 
easier to obtain. The equivalence between time and space averaging is not an 
obvious assumption, but since jams are moving, space and time correlations 
are combined in some way [7] and we consider this assumption to be quite 
safe. In what follows, we will therefore compute the FD by considering either 
the joint probability measure P\(d, 4>) for an open system, or the conditional 
probability measure P\(<j>\d) for a closed system, where 

L number of queues 

L 

N = Xi number of vehicles 

i=i 

L 

<P = ^ i?jl/x 4 >o} integrated flow 

i=l 

are spatial averaged quantities and represent respectively the density and the 
traffic flow. We perform the analysis in the ring geometry: this avoids edge 
effects, fixes the numbers iV of vehicles and L of queues, and finally makes 
sense as an experimental setting. In the statistical physics parlance, the fact 
that N is fixed means that we are working with the canonical ensemble. As a 
result this constraint yields the following form of the joint probability measure: 

ns) = ^ m i{Px{ Xi ,» i ), 

i—l 

with the canonical partition function 

Z L [N] d = f f[Px(xi)S(N - Ef=i Xi). 



d = 



N 

N + L' 
$ 

N + L' 
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A queueing theory approach for a multi-speed exclusion process. 



7 



These expressions are actually independent of A in this specific ring geometry. 
The density-flow conditional probability distribution takes the form 

p ,,,n Z L [N,$] 

with 
and 

Z L {N,<P}^J2 [■[ II i ^( a: *' d wM(^-Ef=ia:0*(*-Ef=iWl{x < >o})- 

{xi} i=l 

(4) 

Note (by simple inspection, see e.g. [5]) that P(<j>\d) is independent of A. 



4 Phase transition and condensation mechanism 

The connection between spontaneous formation of jams and the Bose-Einstein 
condensation has been analyzed in some specific models, with e.g. quenched 
disorder [9], where particles are distinguishable with different but fixed hop- 
ping rates attached to them. In the present situation, all particles are identical, 
but hopping rates may fluctuate, which is related to annealed disorder in sta- 
tistical physics. The condensation mechanism for zero range processes within 
the canonical ensemble has been clarified in some recent work [10]. Let us 
translate in our settings the main features of the condensation mechanism. 
Assume that the number of clients X of an isolated queue has a long-tailed 
distribution 

P(X = x) oc — , a > 1. 
The empirical mean queue size reads 



X = j J2 X u and EX = E X (X) = 



where E\(X) is the expected number of clients in an isolated queue, when 
the arrival rate is A. Within the canonical ensemble, X is fixed, while for the 
grand canonical ensemble, only the expectation E(X) is fixed. In both cases, 
for a > 2 there exists X c such that, when X > X c (E(X) > X c in the grand 
canonical ensemble), one of the queues condenses, i.e. carries a macroscopic 
number of particles. When a > 2, there is a condensate with probability 
weight OIL 1 - 01 ). 

This condensation corresponds to a second order phase transition, and 
occurs at a critical density d c which is the same in the canonical and grand- 
canonical formalism. To determine d c , first consider 
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J(A) d = f EaW 



E\(X) increases monotonically with A, which cannot exceed /io (see Section 2). 
Therefore, if E /10 (X) = x c < oo, then there exists a critical density 

U r . = — , 



such that for d > d c one of the queues condenses. The interpretation is that 
iV c = Lx c is the maximal number of clients that can be in the queues in the 
fluid regime, and the less costly way to absorb the excess N — N c is to put it in 
one single queue. Let us give an example, by specifying the joint law through 

P(»o<R<»oy) = F( fM) y) = (^j) Q , i<y<r, (5) 

where r > 1 is ratio between the highest and lowest speed. In that case, using 
(2) we have the following asymptotic as £ — »■ oo 



P X (X = x) a 



and E^ (X) < oo when a > 2, which yields the possibility of condensation 
above the critical density 

(r - I)"" 2 

^)= Q -2 + (r - I).-* " (6) 



5 Numerical results 

The analysis of (3) in the ring geometry can in principle be performed by 
means of saddle point techniques [10, 11], which we postpone to another work. 
Instead we present a numerical approach: the fd presented in Fig. 5(a)-(c) is 
obtained by solving the recursive relation 

N 

Z L {N,<P] =P X (X = Q)Z L _ 1 {N,<P] + J2 / Px(x,dn)Z L ^[N-x,<P-fit {x>0} }, 

up to some value Lmax = 100 for the number of queues, with a fixed value 
of A < n . Although in principle one arbitrary value of A should suffice, 
in practice, the results for different values of A have to be superposed in 
the diagram to get significant results. Since this recursion is only tractable 
with a finite number of possible velocities, the distribution F used here is 
concentrated to two values fi and \x\. The presence of a discontinuity in the 
fundamental diagram for sm; ill values of r\ = P(R = /j ) /P{R = Hi) is a finite 
size effect, which disappears when the system size is increased while -q is kept 
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Fig. 5. Probabilistic fundamental diagram on the ring geometry (L + N = 100) for 
a two-level speed distribution with rj = 1 (a), 0.001 (b) and 10 -6 (c); corresponding 
single queue distribution as a function of the density for r\ = 0.001 (d). 




Fig. 6. Mean flow as a function of density for a continuous speed distribution (a = 3) 
on the ring geometry with varying sizes L + N (left) and corresponding standard 
deviation rescaled by vL (right). 

fixed. Nevertheless, the direct simulation of the closed L-stage tandem queues, 
with continuous distribution (5), indicates as expected a second order phase 
transition when a > 2 (Fig. 6). This transition is related to the formation 
of a condensate, which is marked by the apparition of a bump in the single 
queue distribution at the critical density (see Fig. 5(d)). This condensation 
mechanism is responsible for the slope discontinuity. Fluctuations scale like 
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l/\/L, as expected from the Central Limit Theorem. Note however that the 
critical density is different than the one given by (6) for the open system. 



6 Perspectives 

In this work, we analyze the fluctuations in the fundamental diagram of traffic 
by considering models from statistical physics and using probabilistic tools. 
We propose a generalization of the tasep by considering a multi-speed ex- 
clusion process which is conveniently mapped onto an L-stage tandem queue. 
When the individual queues are reversible, general results from queueing net- 
work theory let us obtain the exact form of the steady state distribution. This 
measure is conveniently shaped to compute the FD. Depending on the speed 
distribution, it may present two phases, the free-flow and the congested ones, 
separated by a second order phase transition. This transition is associated to 
a condensation mechanism, when slow clusters are sufficiently rare. 

In practice, it is conjectured [8] that there are three phases in the FD, 
separated by first order phase transition. A large number of possible extensions 
of our model are possible, by playing with the definition of the state graph 
of a single queue (Fig. 4(a)). This graph accounts either for the dynamics of 
single vehicle clusters, when queues are associated to empty sites, or to the 
behavior of single drivers when queues are associated to occupied sites. In 
order to obtain first order phase transitions, we will consider in future work 
models where the single queues are not reversible in isolation, for example 
because of an hysteresis phenomenon (Fig. 4(b)). 
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